The effect of self-afRne fractal roughness of wires on atom chips 



Z. Moktadir0 B. DarquieQ M. Kraft0 and E.A. HindsQ 

Atom chips use current flowing in lithographically patterned wires to produce microscopic mag- 
netic traps for atoms. The density distribution of a trapped cold atom cloud reveals disorder in 
the trapping potential, which results from meandering current flow in the wire. Roughness in the 
edges of the wire is usually the main cause of this behaviour. Here, we point out that the edges of 
microfabricated wires normally exhibit self-aflrne roughness. We investigate the consequences of this 
for disorder in atom traps. In particular, we consider how closely the trap can approach the wire 
when there is a maximum allowable strength of the disorder. We comment on the role of roughness 
in future atom-surface interaction experiments. 
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I. INTRODUCTION 

Atom chips are microfabricated structures that allow 
the preparation and manipulation of cold atom clouds or 
Bose-Einstein condensates (BEC) above a substrate sur- 
face. Often, these structures use current-carrying wires 
to produce tightly confining magnetic microtraps close to 
the substrate surface, where atom clouds can be held still, 
transported, or split [l H, 0, Q . With the integration 
of optical components a] and movable structures [3] into 
atom chips, new possibilities are now opening for neutral 
atoms on a chip to form quantum sensors, clocks and in- 
formation processors [8j. However, the homogeneity and 
stability of atom clouds can be compromised close to a 
metallic surface by physical factors that cause fragmen- 
tation and/or the loss of atoms. Two main phenomena 
have been identified: (i) spatial imperfections of the wire, 
which cause the current to flow non-uniformly and make 
the atom trap rough, and (ii) thermal fluctuations of the 
magnetic field near the surface, which drive spin flips of 
the atoms and cause loss [9(. The first of these is the 
subject of our paper. 

Recently, the corrugation of magnetic fields close to 
a wire has been studied extensively. Initial experiments 
showed that atom clouds break up into fragments as they 
approach the surface [13, El , then it was demonstrated 
that this is due to a magnetic field component parallel 
to the wire [l2| , caused by transverse components of the 
current density. Some theoretical efforts were made to 
relate this to the details of the current flow and to 
roughness of the surface and irregularity in the edges of 
the wire [HI, which cause the current to meander. 
If the meander has a single spatial Fourier component of 
wavevector qo, the decay of this anomalous field decreases 
with distance d above the wire according to the Bessel 
function Ki(dqo) to a good approximation When 
the transverse current has a broad noise spectrum, this 
decay can sometimes be described by a power law [12j . 

Conductors lithographically patterned on an atom chip 



are usually fabricated with good bulk homogeneity in 
order to minimise this potential problem of magnetic 
roughness. The width of the wire is typically compa- 
rable to the distance between the magnetic trap and the 
surface. In these typical cases, the meandering of the 
current is driven mainly by the roughness of the edges, 
as determined by the fabrication process. Three pro- 
cesses are available, namely, (i) electrodeposition of the 
metal into a mould formed by a thick photoresist [H, [l5| , 
(ii) etching of a complete metallic film to create the 
space between wires |15j | using wet chemicals or ion beam 
milling [IfJ , and (iii) evaporation of the metal onto the 
substrate through a mask formed by a patterned resist 
(the method known as lift-off). 

With all these methods of microfabrication, the edg es 
of the wires exhibit self-afhne fractal roughness [13, llSj . a 
type of roughness that we now discuss. Consider an edge 
along the z direction with roughness fluctuations f(z). 
By definition the height-height correlation function is 
given by G(r) 2 = ([/(z) — /(z + r)] 2 ), the autocorrelation 
function is C(r) = (f(z)f(z + r)), and the mean square 
roughness is a 2 = (f(z) 2 ). The angle brackets denote av- 
eraging over the (large) length of the wire. These quan- 
tities are connected by the relation G{r) 2 — 2er 2 — 2C(r). 
A self-afRne fractal edge is one that satisfies the scaling 
law G(r) oc r Q , where a is known as the roughness ex- 
ponent or Hurst exponent. The statistical properties of 
such an edge are invariant when the length is scaled by 
a factor A, provided there is an accompanying scaling of 
the transverse dimension by A". Microfabricated edges 
exhibit precisely this type of behaviour on small length 
scales [19j up to a characteristic length £, known as the 
correlation length. For r > £, G{r) tends to the constant 
value y/2a and the autocorrelation function C(r) tends 
to zero. This behaviour is captured by the empirical au- 
tocorrelation function 
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which has the required assymptotic behaviour at large 
and small r and fits experimental data well 20]. The 
Hurst exponent a is normally between and 1 |2l| , while 
the correlation length and rms r oug hness are both typi- 
cally in the range 1-100 nm [l3l 1 1 9L [2fj|. 

Considerable progress has been made in understand- 
ing how roughness in the edges of a wire can gener- 




FIG. 1: Sketch of the ideal wire geometry under consideration. 
Current I flows uniformly along z. At height d above the wire, 
this makes a field along y. Roughness in the edges of the wire 
cause the current to deviate from side to side, producing a 
noise field 8B Z . 



ate r oug hness in the magnetic traps produced by the 
wire [llllT3l.[T3|. However, the analyses to date have con- 
sidered edges with a white noise spectrum or other rather 
specific model spectra. Here we reconsider the magnetic 
noise of atom traps, taking into account this realistic and 
more generally applicable model of the self-affine fractal 
edge roughness. 



II. DESCRIPTION OF SELF-AFFINE 
ROUGHNESS ALONG AN EDGE 

Figure Q] defines a coordinate system and shows the 
wire that we are considering. The left and right edges 
of the wire lie on y = and y = yo with fluctuations 
Jl{z) and Jr(z) respectively. Hence the centre of the 
wire lies on y /2+6y(z), where dy(z) = \ [f L (z) +/«(»], 
with correlation function C(r) = (Sy(z + r)Sy(z)). Wc 
take as our starting point the empirical correlation func- 
tion C(r) given in equation ([1]). In the particular case 
when a = 1/2, the corresponding power spectrum has 
the Lorentzian form 



C{r)e- iqr dr = a 2 £ 



2/tt 



i + q 2 e ' 



(2) 



where q spans the range to oo. In order to have analyt- 
ical results for a more general range of possibilities, we 
extend the power spectrum of equation ([2]) to the form 



P{a,q) = a 2 i 



2/tt 



(1 + aq 2 £ 2 )2 



= a 2 £,P. 



(3) 



This is a one-dimensional version of the approximation 
introduced by Palasantzas (section IV of [20() to describe 
surface noise. The parameter a in the denominator of this 
spectrum is needed to ensure that the integral of equation 
([3|) over all q yields the mean square roughness a 2 . This 
normalisation condition requires 
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FIG. 2: Spectrum of wire roughness P defined by equation 
([3]) versus q(,, where £ is the correlation length. Spectra are 
plotted for two Hurst exponents, a — 1/4 and a = 1. Inset are 
two representative plots showing the roughness over a length 
2£ along z with these two values of a. 



where T is the Euler function. We find by direct numer- 
ical integration that the power spectrum corresponding 
to equation (fTJ) is reasonably well reproduced by equa- 
tion ©, but only over the range 1/4 < a < 1. Figure [5] 
shows the dimensionless spectrum P, defined in equation 
Q, for these two extremes of the Hurst exponent a. At 
low frequency it has the value (2/-7r)(cr 2 £) regardless of a, 
but as the frequency increases, the spectrum with lower a 
also has lower noise. At higher frequency still, this nec- 
essarily reverses because these spectra are normalised. 
Typical profiles of the centre position y /2 + Sy(z), plot- 
ted over a length 2£, are inset into the figure to illustrate 
this. The case of a — 1/4 exhibits more high-frequency 
noise but less long- wavelength noise than that of a = 1 . 

This model is expected to describe fluctuations in the 
centre of a wire fabricated on an atom chip. In par- 
ticular, it describes the atom chip wires currently being 
used at Imperial College London, which are made by ion 
beam milling a gold film. These typically have roughness 
a ~ 3 nm, correlation length £ ~ 20 nm and Hurst expo- 
nent a ~ 0.5. The noise reported in figure 8 of Schumm 
et al. [131] is also consistent with this model, giving for 
their evaporated wire the values a ~ 1.4 nm, a ~ 0.8 
and £ ~ 50 nm. (Note, however that such analysis of 
the power spectrum is not a very reliable way to measure 
£ or oi [22]). In the following, we discuss the roughness 
of the magnetic atom traps produced when current flows 
through such a wire and we investigate how the field fluc- 
tuations vary with the Hurst exponent a and correlation 
length £. Wires made by electrodeposition into a thick 
photoresist mould have so far been much rougher. For 
example, the electroplated wire of [l3[ had a ~ 70 nm, 
a ~ 0.5 and £ ~ 200 nm. Moreover, the spectrum of that 
wire exhibited a second power-law region with exponent 
—2.2 at wavelengths longer than 20 /im, indicating a sec- 
ond regime of correlated roughness. 
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Dimensionless spatial frequency q d 

FIG. 3: Scaling factor f 2 as a function of the dimensionless 
frequency qd. This quantity links the noise in the edges of 
the wire to the noise in the magnetic atom-trapping potential. 
Curves are shown for three ratios of height to width: d/yo 
1 (solid curve), 2 (dashed) and 0.6 (dotted). 




Dimensionless spatial frequency q d 

FIG. 4: Plot of the product Pf 2 , which determines the rough- 
ness spectrum S of the magnetic atom trap. The Hurst expo- 
nent is a — 1 and d = 2yo- Curves are shown for four ratios 
of correlation length £ to distance d: £/d = 33 (dotted line), 
10 (dash-dotted line), 1 (solid line), 0.01 (dashed line). 



III. ROUGHNESS OF THE ATOM TRAP 
FORMED BY A WIRE 

Consider the wire in figure [TJ If the current / flows 
uniformly along the z-direction, the magnetic field lines 
lie in the xy plane. In particular, the field points in the 
y-direction above the centre of the wire. This field is 
cancelled at a height d by applying an opposite uniform 
bias field and the resulting line of zero magnetic field is 
surrounded by a transverse quadrupole field. Magnetic 
atoms can then be trapped at height d above the centre 
of the wire by the magnetic dipole interaction — p, ■ B. 
A small uniform bias field B z is often applied along the 
^-direction as well, so that the magnetic field minimum 
goes to B z rather than to zero. This suppresses the loss 
of atoms through non-adiabatic spin flips. 

In reality, the noise in the edges of the wire causes the 
current to deviate from side to side, generating a noise 
field component 5B Z along the z direction. Consequently, 
the potential energy along the centreline of the trap is no 
longer the constant —fi z B z but is modulated by noise 
~/i z SB z . Assuming that the wire is thin (along x) com- 
pared with the height d, the power spectrum of this noise 
can be written as 

S(q) = f z B 2 P~f = £Bl f^p) S, (5) 

where B = /.i I/2ird is characteristic of the ideal field 
produced by the wire and <J 2 ^P is the power spectrum 
describing fluctuations in the centre of the wire, which we 
take here to be given by equation © . The dimensionless 
scaling factor f 2 translates the noise in the centre of the 
wire to the noise in the field. It is given by [3] 

- = 2sinh(l gyo ) y, 

qy Q sinh(q?/o) ^ n\(2qd) n 

*[i2n+i{jqyo) - i2n+i{-\qyo)}, (6) 



where K n (x) is the modified Bessel function of the second 
kind and r ) n {x) is the incomplete Gamma function. This 
expansion is useful in the range d > yo/2, where a small 
number of terms is sufficient to achieve convergence: 50 
terms at d = 0.6yo and fewer terms at larger distance. 
When d < yo/2, the individual terms become excessively 
large and the series appears not to converge. In summary, 
the spectrum S(q) of the noise in the magnetic atom trap 
depends on the Hurst exponent a and four length scales: 
a and £ in P(a,q), which characterise the roughness of 
the edges and yo and d in /, which define the geometry 
of the trap. The over-all energy scale is given by [i z B§. 

The frequency dependence of f 2 is illustrated in figure 
[3] The three curves correspond to wire widths of d/10 
(solid line), d/2 (dashed) and d/0.6 (dotted). When y is 
small compared with d, f is quite insensitive to its value, 
but / becomes small for a wider wire as the edges move 
further away compared with d. At low frequencies, i.e. 
when qd « 1, the function f 2 ^ [qd (^) arctan(Jg)] 2 

increases in proportion to q 2 , reaching a maximum in 
the vicinity of qd — 1. For large qd, the function de- 
cays as f qd( 2 ^) 2 exp{-qd{2+ (|^) 2 ]}. If the current 
wanders periodically from side to side with a given am- 
plitude, the angular variation of the current density j is 
inversely proportional to the wavelength. Consequently, 
the transverse component j y is proportional to the fre- 
quency q. This is the physical cause of the linear cutoff 
in / at low-frequency. The exponential (Bessel) cutoff 
at high frequency is due to Laplace's equation for the 
magnetostatic potential, which naturally smoothes high 
frequency ripples as one goes far away from the wire. 

The dimensionless spectrum S of the noise in the trap- 
ping potential (see equation ((5|)) is the product of the two 
spectra P and / 2 , shown in figures and O Whereas P 
depends on the frequency through q^, f is a function of 
qd, therefore the shape of the spectrum S depends on 
the ratio d/£ as illustrated in figure 21 The dashed curve 
in figure [?] represents the case of short correlation length, 
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FIG. 5: Roughness spectra S of the magnetic trap for two 
values of the Hurst exponent a. Light curves: a = 1. Heavy 
curves: a = 1/4. As in figureUl spectra are given for £/d = 33 
(dotted), 10 (dash-dotted), 1 (solid), and 0.01 (dashed). 

£ = d/100, for which g£ << 1 over the whole range of the 
graph, making P constant at = 2 /it. In this limit, equa- 
tion (|5|) gives S(q) — (J-IBq ( ^jfj f 2 , a spectrum that is 

independent of the wire roughness except for the cr 2 £ in 
the over-all scale factor. This general behaviour of a spec- 
trum proportional to f 2 persists throughout the range 
£ < d, as also indicated by the solid line representing 
£ = d. At the other extreme, the dotted curve in figure 
2] represents the case of long correlation length, £ = 33d, 
for which qd <^ 1 over the whole range of interest, giv- 
ing f 2 ~ (qd) 2 and therefore S(q) ~ fi 2 B 2 a 2 ^q 2 P. This 
general behaviour of S(q) oc q 2 P is characteristic of the 
whole range £ > d, as also illustrated by the dash-dotted 
line in figure [4j 

Figure [5] shows the effect of changing the Hurst expo- 
nent in the roughness spectrum of the wire from a — 1 , 
as in figure d] (light curves), to a = 1/4 (heavy curves). 
When £ <C d (dashed curves), the change of Hurst expo- 
nent makes no difference because the spectrum is essen- 
tially independent of P. By contrast, the dotted curves 
representing £ = 33c£ exhibit a strong dependence on the 
Hurst exponent. Reducing a from 1 to 1/4 suppresses the 
low frequency noise and increases the power at higher fre- 
quencies, as already noted in the context of figure[21 This 
has the effect of moving the peak of the noise spectrum to 
higher frequencies. The same effect is seen in the dash- 
dotted curves of figure [5] representing £ = Wd. When £ 
is equal to d (solid curves), the change to a = 1/4 sup- 
presses the low frequency part of the spectrum, but the 
corresponding increase at higher frequency is not evident 
because the spectrum is cut off at higher frequencies by 
the exponential roll-off of the function f 2 . 

The mean square roughness of B z , let us call it V for 
variance, is obtained by integrating equation ([5]) over fre- 
quency: 

1 f°° 

V(a, e, d, y ) = (B 2 ) = — / S(q)dq . (7) 
A* 2 Jo 

Once again, it is useful to adopt a dimensionless version 
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FIG. 6: Dimensionless magnetic field variance V versus dis- 
tance d of the trap from the wire, normalised to the width 
j/o of the wire. Curves are given for f/d = 1 (solid), 20 (dot- 
ted) and 0.01 (dashed). Light curves: a = 1. Heavy curves: 
a = 0.25. 

in order to understand how V depends on the various 
length scales involved. The form 

n£,M = (~) Wo =l SZdq (8) 

renders the field variance dimensionless and makes it a 
function of d/yo, d/£ and a. This is plotted versus d/yo 
in figure [S] for three values of d/t; and two values of a. 
We see immediately that the dimensionless variance V 
is approximately constant with distance, corresponding 
to a l/d A variation in V that weakens when d < yo. 
This result is consistent with experimental observations 
that the noise decreases with increasing distance from 
the wire [ill, EH EH • All the dependence on £ and a is 
contained in V, which is largest when d ~ £ (solid lines). 
Changing a from a = 1 (light curves) to a = 1/4 (heavy 
curves) makes essentially no difference when £/d is small 
(dashed lines), because P is effectively a constant under 
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FIG. 7: Dimensionless magnetic field variance V versus the 
ratio d/£ for the particular case d = yo- Solid line: a = 1, 
dashed line: a = 1/2, dotted line: a = 1/4. Dashed-dotted 
line: f/d law given by equation © and valid when d>^. 
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these conditions, as discussed above. By contrast, the 
same change of a reduces the noise when d ~ £ (solid 
line) and increases it when d <C £ (dotted line). Another 
view of the same parameter space is given in figure [7J 
which plots the field noise V at height d = yo versus 
d/£, for three values of a. This shows more clearly the 
peaking of the noise near £ = d, and the insensitivity to 
the value of the Hurst exponent when £ is small. 

The best present methods for fabricating atom chip 
wires yield correlation lengths of order 1-100 nm, while 
current experiments operate at distances in the range 1— 
100 pxa. This places experiments firmly in the domain 
of small £/d, where the value of the Hurst exponent does 
not influence the roughness of the magnetic trap signifi- 
cantly and we can take P ~ 2/n. In this regime we find 
that V oc £ I d. For any particular value of yo we can in- 
tegrate equation ([8J) numerically to obtain the constant 
of proportionality. For example, with y — d we find 



l,f » l,a) 



0.274 i . 
d 



(9) 



For narrower wires, i.e. for d > yo, this constant changes 
very little as we have already seen in figure [6] 



IV. CONSEQUENCES FOR THE DESIGN OF 
THE ATOM TRAP 



One of the primary motivations for atom chips is to 
achieve small traps with high trapping frequencies. This 
requires a high field gradient, which is approximately 
Hoi /(2ird 2 ) for narrow wires (by which we mean d < yo). 
The tightest confinement is achieved by bringing the 
atoms close to narrow wires, even though smaller wires 
carry less current [l3j , but this also increases the rough- 
ness of the magnetic field. In many experiments, the field 
has a maximum permissible variance, let us call it V max . 
For example, a Bose-Einstein condensate will break into 
separate clouds unless the noise in the trapping potential 
is smaller than the chemical potential. In these cases, 
the limit on trap roughness imposes a minimum distance 
from the surface and hence a maximum achievable mag- 
netic field gradient. 

The heat generated by electrical resistance limits the 
current that can be tolerated in a lithographically fab- 
ricated wire to I max = ^yoy/xo [HI, where xo is the 
thickness of the wire (figure [1]) and the constant K char- 
acterises the heat flow across the interface between the 
wire and the substrate. When the maximum current is 
passed through a wire having yo — d, equations JS]) and 
© give the variance of B z as 



V = 0.274 



d 3 



2vr 



xo 



(10) 



of closest approach and the maximum field gradient are 

1/3 



d rtl 



0.274 a^(^ 

Holm ax 



\ 2n J V ma . 

HoKy/XoVm 



2nd: 



0.274 x 2ncr 2 Z 



1/3 



(11) 
(12) 



Let us take \/V max = 1 mG, for which the rms roughness 
in the potential corresponds to a temperature of 67 nK. 
We take k ~ 3 x 10 7 A m~ 3 / 2 0, Sg|, which is typical 
for gold wires on a Si/Si0 2 substrate. We further assume 
the values xq ~ 1 /im, a ~ 3 nm, and £ ~ 20 nm, which 
are typical of our present wires, as discussed in sectionHTl 
Then the minimum distance of approach to the wire is 
d m i n — 6 /im, the current in the wire is I max — 170 mA, 
and the corresponding maximum field gradient is B' max ~ 
11 Tern -1 . 

In the presence of a bias field B z along the z-direction 
(figure [TJ, the potential energy of the trapped atom near 
its equilibrium position is Hz(B 2 + (B^^p) 2 ) 1 ^ 2 , where p 
is the transverse (xy) displacement. The corresponding 
frequency for transverse harmonic oscillations is 



fn 



1 

2tt~ 



B' 



mB z 



(13) 



Upon setting this equal to V max , we find that the distance 



where m is the mass of the atom. Taking a typical 
value of B z ~ 0.5 G, the maximum transverse frequency 
for 87 Rb atoms in the F = 2, mp = +2 ground state 
is fmax — 190 kHz. This result indicates that atom 
chips can achieve very high trapping frequencies, com- 
parable with those already demonstrated in optical lat- 
tices, while remaining adequately smooth. The cor- 
respondingly small extension of the vibrational ground 
state wavepacket is only, 17 nm, making such traps very 
promi sing for studing the physics of 1-dimensional cold 
gases [241 ] - 



V. DISCUSSION 

We have examined how a magnetic trap formed by a 
current-carrying wire is sensitive to the roughness on the 
edges of the wire. In particular, we have extended previ- 
ous discussions to consider the case of a self-affine fractal 
roughness spectrum with a correlation length £. This is 
of interest because the methods used to fabricate wires 
on an atom chip generally produce such roughness. Our 
analysis has shown how the spectrum of trap roughness 
involves an interplay between the spectrum of the wire 
roughness and the spectrum of the transfer function that 
converts deviations of the current into fluctuations of the 
magnetic field. In most current experiments, there is a 
clear hierarchy of length scales in which d ~ yo 3> £. 
When there is also a maximum acceptable roughness of 
the magnetic trap, this leads to equation {TTJ for the min- 
imum operating distance between the trap and the wire. 
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There is correspondingly a maximum achievable field gra- 
dient given by equation (Tr2"|) . These results argue for 
minimising the quantity cx 2 £ because this determines the 
spectral density of the edge roughness at low frequency, 
which is what generates the noise in the magnetic trap. 
In these cases where d ~ i/q ^> £, the Hurst exponent 
is less important because it only affects the spectrum at 
frequencies above l/£, which do not contribute signif- 
icantly to the trap roughness. Naturally, if there is a 
second power- law regime at long wavelength, as in some 
wires [131 ] , then the corresponding Hurst exponent signif- 
icantly affects the magnetic trap roughness at that length 
scale. 

There is considerable interest in trapping atoms much 
closer to the surface. Atoms trapped at sub-micron dis- 
tances would begin to probe the details of the short-scale 
noise and in that case the value of a would be significant. 
In the range of 0.1-1 fim one could intentionally create 
surfaces with a variety of noise spectra in order to prop- 
agate BEC through custom-made disorder potentials as 
a study of quantum localisation phenomena, as proposed 
by [IH ■ Experiments of this type have already been done 
using random optical potentials [H|,[26|,[23], but this mag- 
netic disorder offers a different noise spectrum with the 
possibility of very short correlation length. The effective 
amplitude of the noise can be conveniently controlled by 
appropriate modulations of the currents that form the 
trap, as recently demonstrated by [28|. A natural lower 
limit to the distance of closest approach is set by the elec- 
tromagnetic attraction of the atom towards the surface 
- the Van der Waals (Casimir-Polder) force — which 
grows as 1/d (1/d 5 ) and overwhelms the trapping force 
at a distance of order 100 nm HI,!!!. This force itself is 
also of fundamental interest and can be measured further 
away from the surface by means of cold atoms [10, Hl| . 
Even more exotically, one can hope to measure the grav- 
itational attraction at short range, which might exhibit 
departures from the Newtonian law as a result of extra 



dimensions [32]. These experiments also require careful 
control over the noise of the surface. We note that if the 
distance to the surface becomes much less than the width 
of the wire, the corrugation of the surface [l3T ] and im- 
perfections of the bulk [33| may contribute significantly 
to the magnetic trap roughness. 

In this study we have not considered the role of ad- 
ditional technical noise at long correlation lengths, but 
such noise certainly exists as a result of imperfections in 
the fabrication process and can make a significant con- 
tribution to the roughness of the magnetic trap at long 
wavelength. For example, a supposedly straight line may 
be bent by optical aberrations during the lithographic 
patterning. A bend of only 100 /irad amounts to a devia- 
tion of only 10 nm over a length of 100 fj,m. Nevertheless, 
in a bias field of Bo = 10 G this would generate an appre- 
ciable unwanted B z of 1 mG. There can also be a periodic 
wobble of straight lines due to the imperfection of me- 
chanical translation stages. These kinds of defects are 
very hard to measure by standard microscopy because 
they involve small transverse displacements over length 
scales larger than the normal field of view in an SEM or 
AFM microscope. In fact, this noise is best measured by 
the cold atoms themselves through its effect on their den- 
sity distribution in the trap. Ref. [33[ has already noted 
that on certain length scales, this method may provide a 
uniquely sensitive way to probe the magnetic field near 
surfaces. 
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